Simulation of a Single Polymer Chain in Solution by Combining Lattice Boltzmann and Molecular Dynamics

In this paper we establish a new efficient method for simulating polymer-solvent systems which combines a lattice Boltzmann approach for the fluid with a continuum molecular dynamics (MD) model for the polymer chain. The two parts are coupled by a simple dissipative force while the system is driven by stochastic forces added to both the fluid and the polymer. Extensive tests of the new method for the case of a single polymer chain in a solvent are performed. The dynamic and static scaling properties predicted by analytical theory are validated. In this context, the influence of the finite size of the simulation box is discussed. While usually the finite size corrections scale as L^{-1} (L denoting the linear dimension of the box), the decay rate of the Rouse modes is only subject to an L^{-3} finite size effect. Furthermore, the mapping to an existing MD simulation of the same system is done so that all physical input values for the new method can be derived from pure MD simulation. Both methods can thus be compared quantitatively, showing that the new method allows for much larger time steps. Comparison of the results for both methods indicates systematic deviations due to non-perfect match of the static chain conformations.Comment: 17 pages, 12 figures, submitted to J. Chem. Phy

Screening of Hydrodynamic Interactions in Semidilute Polymer Solutions: A Computer Simulation Study

We study single-chain motion in semidilute solutions of polymers of length N = 1000 with excluded-volume and hydrodynamic interactions by a novel algorithm. The crossover length of the transition from Zimm (short lengths and times) to Rouse dynamics (larger scales) is proportional to the static screening length. The crossover time is the corresponding Zimm time. Our data indicate Zimm behavior at large lengths but short times. There is no hydrodynamic screening until the chains feel constraints, after which they resist the flow: "Incomplete screening" occurs in the time domain.Comment: 3 figure

A new model for simulating colloidal dynamics

We present a new hybrid lattice-Boltzmann and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of spherical colloidal particles. The solvent is modeled on the level of the lattice-Boltzmann method while the molecular dynamics is done for the solute. The coupling between the two is implemented through a frictional force acting both on the solvent and on the solute, which depends on the relative velocity. A spherical colloidal particle is represented by interaction sites at its surface. We demonstrate that this scheme quantitatively reproduces the translational and rotational diffusion of a neutral spherical particle in a liquid and show preliminary results for a charged spherical particle. We argue that this method is especially advantageous in the case of charged colloids.Comment: For a movie click on the link below Fig

A Numerical Model for Brownian Particles Fluctuating in Incompressible Fluids

We present a numerical method that consistently implements thermal fluctuations and hydrodynamic interactions to the motion of Brownian particles dispersed in incompressible host fluids. In this method, the thermal fluctuations are introduced as random forces acting on the Brownian particles. The hydrodynamic interactions are introduced by directly resolving the fluid motions with the particle motion as a boundary condition to be satisfied. The validity of the method has been examined carefully by comparing the present numerical results with the fluctuation-dissipation theorem whose analytical form is known for dispersions of a single spherical particle. Simulations are then performed for more complicated systems, such as a dispersion composed of many spherical particles and a single polymeric chain in a solvent.Comment: 6 pages, 8 figure

Alkali and Alkaline Earth Metal Compounds: Core-Valence Basis Sets and Importance of Subvalence Correlation

Core-valence basis sets for the alkali and alkaline earth metals Li, Be, Na, Mg, K, and Ca are proposed. The basis sets are validated by calculating spectroscopic constants of a variety of diatomic molecules involving these elements. Neglect of $(3s,3p)$ correlation in K and Ca compounds will lead to erratic results at best, and chemically nonsensical ones if chalcogens or halogens are present. The addition of low-exponent $p$ functions to the K and Ca basis sets is essential for smooth convergence of molecular properties. Inclusion of inner-shell correlation is important for accurate spectroscopic constants and binding energies of all the compounds. In basis set extrapolation/convergence calculations, the explicit inclusion of alkali and alkaline earth metal subvalence correlation at all steps is essential for K and Ca, strongly recommended for Na, and optional for Li and Mg, while in Be compounds, an additive treatment in a separate `core correlation' step is probably sufficient. Consideration of $(1s)$ inner-shell correlation energy in first-row elements requires inclusion of $(2s,2p)$ `deep core' correlation energy in K and Ca for consistency. The latter requires special CCV$n$Z `deep core correlation' basis sets. For compounds involving Ca bound to electronegative elements, additional $d$ functions in the basis set are strongly recommended. For optimal basis set convergence in such cases, we suggest the sequence CV(D+3d)Z, CV(T+2d)Z, CV(Q+$d$)Z, and CV5Z on calcium.Comment: Molecular Physics, in press (W. G. Richards issue); supplementary material (basis sets in G98 and MOLPRO formats) available at http://theochem.weizmann.ac.il/web/papers/group12.htm

Local Simulation Algorithms for Coulomb Interaction

Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte-Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow equilibration using a local Hamiltonian. The method introduces an auxiliary field with constrained dynamics so that the equilibrium distribution is determined by the Coulomb interaction. We demonstrate the efficiency of the method by simulating a simple, charged lattice gas.Comment: Last figure changed to improve demonstration of numerical efficienc

How does the 'ancient' asexual Philodina roseola (Rotifera:Bdelloidea) handle potential UVB-induced mutations?

Like other obligate asexuals, bdelloid rotifers are expected to suffer from degradation of their genomes through processes including the accumulation of deleterious mutations. However, sequence-based analyses in this regard remain inconclusive. Instead of looking for historical footprints of mutations in these ancient asexuals, we directly examined the susceptibility and ability to repair point mutations by the bdelloid Philodina roseola by inducing cyclobutane-pyrimidine dimers (CPDs) via exposure to UVB radiation (280-320 nm). For comparison, we performed analogous experiments with the facultative asexual monogonont rotifer Brachionus rubens. Different strategies were found for the two species. Philodina roseola appeared to shield itself from CPD induction through uncharacterized UV-absorbing compounds and, except for the genome reconstruction that occurs after desiccation, was largely unable to repair UVB-induced damage. By contrast, B. rubens was more susceptible to UVB irradiation, but could repair all induced damage in similar to 2 h. In addition, whereas UV irradiation had a significant negative impact on the reproductive output of P. roseola, and especially so after desiccation, that of B. rubens was unaffected. Although the strategy of P. roseola might suffice under natural conditions where UVB irradiation is less intense, the lack of any immediate CPD repair mechanisms in this species remains perplexing. It remains to be investigated how typical these results are for bdelloids as a group and therefore how reliant these animals are on desiccation-dependent genome repair to correct potential DNA damage given their obligate asexual lifestyle.</p

Algorithm for numerical integration of the rigid-body equations of motion

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.Comment: 4 pages, 1 figur

Dynamics and Scaling of 2D Polymers in a Dilute Solution

The breakdown of dynamical scaling for a dilute polymer solution in 2D has been suggested by Shannon and Choy [Phys. Rev. Lett. {\bf 79}, 1455 (1997)]. However, we show here both numerically and analytically that dynamical scaling holds when the finite-size dependence of the relevant dynamical quantities is properly taken into account. We carry out large-scale simulations in 2D for a polymer chain in a good solvent with full hydrodynamic interactions to verify dynamical scaling. This is achieved by novel mesoscopic simulation techniques

Electron correlations for ground state properties of group IV semiconductors

Valence energies for crystalline C, Si, Ge, and Sn with diamond structure have been determined using an ab-initio approach based on information from cluster calculations. Correlation contributions, in particular, have been evaluated in the coupled electron pair approximation (CEPA), by means of increments obtained for localized bond orbitals and for pairs and triples of such bonds. Combining these results with corresponding Hartree-Fock (HF) data, we recover about 95 % of the experimental cohesive energies. Lattice constants are overestimated at the HF level by about 1.5 %; correlation effects reduce these deviations to values which are within the error bounds of this method. A similar behavior is found for the bulk modulus: the HF values which are significantly too high are reduced by correlation effects to about 97 % of the experimental values.Comment: 22 pages, latex, 2 figure